function [f,P1] = myFFT(input, Fs)
% [f,P1] = myFFT(input, Fs)
% 快速傅里叶变换，输入时间序列，分解为频域
% 输入变量：
% 	input：待处理的信号
%	Fs：	取样频率
% 输出变量：
%	f: 分解信号的频率
%	P1：分解信号的振幅
% 另请参阅：
%	myFFTshow()
%	FFT_tutorial.m

% 作者：sid（cnqdcyq@163.com）
% 日期：2021/7/28
L = length(input);
Y = fft(input);
P2 = abs(Y/L); % 计算双侧频谱
P1 = P2(1:L/2+1); 
P1(2:end-1) = 2*P1(2:end-1); % 计算单侧频谱

f = (Fs * (0:(L/2)) / L)';
end
